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Abstract 

A current can be induced in a closed device by changing control parameters. The amount Q of particles that are transported 
via a path of motion, is characterized by its expectation value (Q), and by its variance Var(Q). We show that quantum 
mechanics invalidates some common conceptions about this statistics. We first consider the process of a double path crossing, 
which is the prototype example for counting statistics in multiple path non-trivial geometry. We find out that contrary to 
the common expectation, this process does not lead to partition noise. Then we analyze a full stirring cycle that consists 
of a sequence of two Landau-Zener crossings. We find out that quite generally counting statistics and occupation statistics 
become unrelated, and that quantum interference affects them in different ways. 
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1. Introduction 

Consider a closed isolated quantum system, say a 
3 site ring as in Fig.l. Quite generally, in the absence 
of magnetic field, the stationary states of the system 
carry zero current. If one wants to have a non-zero cur- 
rent I through a section of the device, one has either to 
prepare it in a non-stationary state or to drive the sys- 
tem. Driving means changing some parameters in time. 
During a time period t the amount of particles that 
get through the section is Q. One may ask what is the 
distribution of the measured Q, and in particular what 
is the expectation value (Q), and what is the variance 
Var(Q). This is known as counting statistics [3,4,5]. 
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Typically the driving is periodic, and Q is defined as 
the amount of particles that are transported per pe- 
riod. The feasibility to have non-zero Q (non zero "DC" 
current) due to periodic ("AC") driving is known in 
the context of open geometry as "quantum pumping" 
[1,2]. We use the term "quantum stirring" [6,7] in or- 
der to describe the analogous effect with regard to a 
closed device [8]. 

The theory of quantum pumping and counting 
statistics in open geometries is well studied. Any at- 
tempt to adopt ideas from this literature to the present 
context of quantum stirring is dangerous, and likely to 
result in major misconceptions, which following [9,10] 
we would like to highlight. For this purpose let us 
consider the simplest model that can be imagined: a 
particle in a 3-site system (Fig.l). We assume that 
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Fig. 1. Toy models that are analyzed in this paper: a par- 
ticle in a 2 site system (upper illustration), and a particle 
in a 3 site system (lower illustration). Initially the parti- 
cle is prepared in the |0) site where it has a potential en- 
ergy u. The hopping amplitudes per unit time (the cs) are 
also indicated. In the case of a 3 site system, the time units 
are chosen such that the hopping amplitude per unit time 
between |1) and 2) equals unity, while the other ampli- 
tudes are assumed to be small (jci|, |c2 1 -C 1). The current 
is measured through the dotted section. 



we have a control over the potential u of the left site, 
and over the couplings ci and C2 that bind the left 
site to the right sites. The transported Q is measured 
through one of these two bonds. Let us introduce two 
simple questions that will be answered later on. 

Question 1. — We start with a very negative u and 
prepare a particle in the left site. Then we gradually 
raise u so as to have an adiabatic transfer of the par- 
ticle to the right side. If we had ci 7^ but C2 = one 
obviously expects in the strict adiabatic limit Q = 1 
with zero variance. We ask: what would be the corre- 
sponding result if we have |ci| = |ca|, and more gener- 
ally how does the result depend on the relative size of 
the couplings? If the |ci| coupling is larger, does the 
result reflect having (say) Q = 1 with 70% probability 
and Q = with 30% probability? 

Question 2. — During a cycle a conventional pump 
takes an electron from the left lead and ejects it to the 
right lead. Hence the pumped charge per cycle for a 
leaky pump is Q < 1. Now we integrate this quantum 
pump into a closed circuit and operate it. We ask what 
is the statistics of Q in the new (integrated) configura- 
tion. Can it be much different? What are the maximal 
(Q) and minimal Var(Q) that can be achieved per cy- 
cle? Are they both proportional to the number of the 
cycles as in the classical reasoning, leading to ~\/t sig- 
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Fig. 2. The adiabatic levels of the 3-site Hamiltonian dur- 
ing one period of a pumping cycle. In the absence of cou- 
pling (ci = C2 = 0) the Eq = u(t) level intersects the sym- 
metric E+ = 1 level. With non zero coupling these inter- 
sections become avoided crossings, and the particle follows 
adiabatically the thickened lines. For presentation purpose 
we indicate that either ci or C2 equal zero ("blocked"), but 
in the general analysis we allow any splitting ratio, includ- 
ing the possibility c\ = C2 of having the same amplitude to 
take either of the two paths. 

nal to noise ratio? 

Outline. — In the next section we define the model 
and the counting operator Q(t). In sections 3 and 4 
we discuss the restricted quantum-classical correspon- 
dence that applies to the analysis of single path cross- 
ing, while in sections 5 we consider multiple path ge- 
ometries. How to treat interference in a sequence of 
Landau-Zener crossings and the analysis of quantum 
stirring are discussed in section 6. The long time count- 
ing statistics is discussed in sections 7 and 8. The rela- 
tion to the theory of spreading and dissipation is illu- 
minated in section 9. The main observations are sum- 
marized in section 10. 



2. Modeling 

We consider the 3 site system of Fig.l which is de- 
scribed by the Hamiltonian 
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(1) 



We label the sites by i — 0, 1, 2. We have control over 
the potential Eo = u of the left site (#0). The two right 
sites (#1 and #2) are permanently coupled to each 
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other forming a double well with energy levels E- — — 1 
and E+ = +1. We also have relatively small couplings 
(|ci| 2 , | C2 1 2 <JC 1) that allow transitions between the left 
and the right sites. 

Later we consider the stirring cycle which is de- 
scribed in Fig. 2. Its operation is inspired by the com- 
mon peristaltic mechanism. Namely, the coupling con- 
stants ci and C2 are regarded as "valves". In the first 
half of the cycle C2 = and u is raised, leading to an adi- 
abatic transfer from-left-to-right via the i— » 1 bond. 
In the second half of the cycle c\ = and u is lowered, 
leading to an adiabatic transfer from-right-to-left via 
the i — * 2 bond. The net effect is to pump one particle 
per cycle. 

The matrix representation of the operator which is 
associated with the current through the i— > 1 bond is 
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In the Heisenberg picture the time dependent current 
operator is defined as T(t) = U(tyIU(t), where U(t) 
is the evolution operator. Consequently the counting 
operator is defined as 



Q = f I{t')dt' 
Jo 



(3) 



The counting operator, unlike the current operator, is 
not a conventional observable in quantum mechanics. 
What can be measured in practice are only the first 
two moments of Q, which are (Q) and Var(Q). Still, 
on the mathematical side, it is convenient to treat Q 
the same way as one treats conventional observables. 
Namely, to regards its eigenvalues Q r as the possible 
outcomes of a measurement, and to associate with a 
given preparation ip a probability distribution 



P(Qr) 
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(4) 



We shall refer to Eq.(4) as the naive definition of the 
full counting statistics. The physical definition of P(Q) 
is much more complicated [3,4], but it leads to the 
same first and second moments [9] . In the present paper 
we are not interested in the full counting statistics, 
but only in the first two moments, and therefore, for 
presentation purpose, we adopt the naive point of view. 



3. Single path crossing 

Let us consider first the very simple case of 2-site 
systems with left site (#0) and right site (#1). The 
particle is prepared in the left site and after some time 
has some probability p to be found in the right site. This 
probability can be regarded as the expectation value 
of the occupation operator N that has the eigenvalues 
or 1 depending on whether the particle is in state #0 
or state #1. 

In the classical analysis the possible outcomes of a 
measurement are N = Q = 1 if the particle goes from 
left to right, and N = Q = otherwise. Obviously one 
should find out that 



P(Q) = ^ 



p 
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for Q = 1 
for Q = 
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It follows that the feth moment is (Q k ) = p and there- 
fore 



(Q) = p 

Var(Q) = (l-p)p 



(6) 
(7) 



In the quantum mechanical treatment the counting 
operator is the integral over the current operator which 
has zero trace, so it should also be a 2x2 traceless 
Hermitian matrix which we write in the i=0, 1 basis as 



Q = 



+Q\\ iQ± 
-iQ* ± -Qii 



(8) 



If the particle is initially prepared in the #0 site then 
(Q) = Q\\ and Var(Q) = |Q±| 2 . So now the question 
arises whether one should expect an agreement be- 
tween the quantum result and the classical results for 
the first and second moments. In the next section we 
argue that for a single path geometry the answer is 
yes: The first and second moments of Q should both 
equal p. It follows that the matrix that represents the 
counting operator in the i=0, 1 basis is expressible in 
terms of p and an extra phase: 



Q = 



+p 



iyj (l-p)p e 10 
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Later we are going to use this expression as a building 
block in the analysis of a multiple crossing scenario. 
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For completeness we note that from the above anal- 
ysis it follows that the eigenvalues of the counting op- 
erator are Q± = ±^/p with 



P(0) 



l (l+v/p)/2 
\ (l-VP)/2 



for Q = Q- 
for Q = Q 4 



(10) 



which should be contrasted with Eq.(5). The proper 
analysis [9] of the full counting statistics gives a more 
complicated quasi-distribution that neither agrees with 
the naive nor with the classical result, but still has the 
same first and second moments. 

The above analysis has assumed nothing about the 
detailed form of the 2x2 Hamiltonian. We merely had 
assumed that the particle was prepared in site #0, and 
that after time t there is some probability p to find the 
particle in site #1. However, in the later sections we are 
going to discuss specifically adiabatic crossings. Dur- 
ing such a process the on-site energies (Eo=u(t) and 
_Bi=const) cross each other, and due to the inter-site 
coupling (=c) the particle is adiabatically transferred 
from #0 to #1. Still there is some small probability 
for a non-adiabatic transition, so called Landau-Zener 
transition [11], that would leave the particle in #0. 
There is a well known formula, that allows to calculate 
the probability of such transition: 



Plz = exp 



- 2 4 

u 



(11) 



Consequently the probability to find the particle in the 
right side at the end of the process is 

V = I -Plz (12) 
Note that it becomes 100% in the strict adiabatic limit. 



4. Quantum-classical correspondence 

For a single path transition we can prove that the 
first two moments of Q should be in agreement with 
the classical expectation. This is based on the relation 
between the occupation operator TV (whose eigenval- 
ues are and 1), and the counting operator Q (whose 
eigenvalues do not have a classical interpretation) . The 
relation is implied by the Heisenberg picture equation 
of motion: 



-Nit) = l(t) 
Integrating over time we get 
Af(t)-Af(0) = Q 



(13) 



(14) 



Assuming that the particle is initially prepared in the 
left site we get for the k — 1, 2 moments of Q 
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(15) 



It is important to realize that the derivation cannot 
be extended to the k > 2 moments because A/"(i) does 
not commute with jV(0). In fact it is not difficult to 
calculate the higher moments: one just has to realize 
that from Eq.(9) it follows that Q 2 = pi, and conse- 
quently the even moments are p k ^ 2 , while the odd mo- 
ments are p( fe+1 )/ 2 . Optionally this result can be ob- 
tained from Eq.(10) as in [9]. 

In the single path transition problem we say that 
we have restricted quantum-classical correspondence. 
The first two moments come out the same as in the 
classical calculation. Encouraged by this observation 
let us speculate what should be the results in more 
complicated circumstances that involve multiple path 
geometries... 

Let us consider our 3 site system (Fig.l). We would 
like to analyze the first half of the cycle which is de- 
scribed by Fig. 2. The particle is initially prepared in 
the left site. We gradually raise u so as to have an adi- 
abatic transfer of the particle to the right side. The oc- 
cupation probability of the right side at the end of the 
process is denoted by p. But now we have to remember 
that the particle could get there either via the i — » 1 
bond or via the 0^2 bond. Motivated by a stochas- 
tic point of view one may argue that the process is like 
partitioning of a current, and therefore 



(Q) = Ap 
Var(Q) = (1 — Ap)Ap [stochastic] 

where the splitting ratio is 



Cl 



Cl' 2 + C2 



[stochastic] 



(16) 
(17) 



(18) 



If for example we have a strict adiabatic process with 
p= 1 and the splitting ratio is A = 1/2, then the 
stochastic expectation is to have Var(Q) = ( 1 /2) 2 . Fur- 
thermore, considering a multi-cycle stirring process, 
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one may argue that the variance should be accumu- 
lated in a stochastic manner: 



Var(Q) oc time 



[stochastic] 



(19) 



All the results above that are labeled as 'stochastic' 
might apply in the case of a non-coherent processes, 
or in the case of open geometries with leads attached 
to equilibrated reservoirs. Below we are going to show 
that the above 'stochastic' results do not apply to the 
analysis of coherent transport in multiple path closed 
geometries. 



5. Double path adiabatic crossing 

For ci = C2 = the 3 site Hamiltonian of Eq.(l) is 
diagonal in the \Eo), \E-}, \E + ) basis. If we limit our- 
selves to processes which involve adiabatic crossings of 
\Eq) and \E+), as in Fig. 2, then transitions to \E-) can 
be neglected and we can work with a 2 x 2 Hamiltonian 
in the \Eo) and \E + ) representation 



n = 




(20) 



where c = (a+C2)/v / 2. In the same representation the 
current operator of Eq.(2) takes the form 



ic 



-ic 



where the splitting ratio is 

Cl 



A = 



Cl + C 2 



(21) 



(22) 



which should be contrasted with the stochastic expres- 
sion Eq.(18). If we had A = 1 it would be the same 
problem as single path crossing. The multiple path ge- 
ometry is reflected in having 1 1— » XT and consequently 
Q I— » XQ. It follows that 



(Q) = Ap 
Var(Q) = A 2 (l-p)p 



(23) 
(24) 



The latter expression for the variance should be con- 
trasted with the stochastic expression Eq.(17). We now 



turn to discuss two surprises that are associated with 
the above results. 

The first surprise comes out from Eq.(22): it is the 
possibility to have A outside of the range [0, 1]. This 
happens if ci and C2 are opposite in sign and close in 
absolute value. So we can cook a cycle such that the 
splitting is into 700% in one path and —600% in the 
other. What does it mean? After some reflection one 
realizes that a proper way to describe the dynamics is 
to say that the driving it induces a circulating current 
in the system. Using a classical-like phrasing one may 
argue that during the transition the particle can encir- 
cle the ring 6 times before it makes the final crossing 
to the right side. 

The second surprise comes out from Eq.(24): it is 
the way in which A appears in the variance calculation. 
Consider for example a strict adiabatic process with 
p = 1. If say A = 1/2 we do not get Var(Q) = (1/2) 2 
but rather Var(Q) = 0. One may say that we do not 
have an incoherent partitioning of the current, but 
rather a noiseless exact splitting of the wavepacket. 
This should be contrasted with the common picture of 
getting shot noise due to partition of current in open 
geometries. 



6. Quantum stirring 

We can use the results that have been obtained in 
previous sections in order to calculate Q for the full 
stirring cycle which is described in Fig. 2. For this pur- 
pose we regard the full stirring cycle as a sequence of 
two Landau Zener crossings, where the first is char- 
acterized by a splitting ratio Ao while the second is 
characterized by a different splitting ratio Ao . The net 
effect is 



<Q> = A -A + O(Plz) 



(25) 



An optional way to derive this result is to make a full 
3 level calculation using the Kubo formula [8] . Here we 
have bypassed the long derivation by making a reduc- 
tion to a multiple path crossing problem. The naive ex- 
pectation is to have \ Q\ < 1 if the valves are leaky. But 
by playing with the splitting ratio we can get \Q\ ^> 1 
per cycle. In the language of the Kubo formalism [8,6] 
this happens if the pumping cycle encloses a degener- 



5 



acy. A large Q reflects a huge circulating current which 
is induced by the driving. 

We can regard the stirring as an induced persistent 
current. Having figured out what is (Q), the next chal- 
lenge in line is to calculate the variance Var(Q). For 
this purpose we regard the full stirring cycle as a se- 
quence of two Landau Zener crossings. The one period 
evolution operator can be written as 



[/(cycle) = [T (7 l ° z T] e"* [U° z ] 



(26) 



We now explain the ingredients of this expression. 
The adiabatic approximation for the I7lz of a single 
Landau-Zener crossing is well known (see e.g. [12]): 



t/ L z 



V1-Plz VK^ 



(27) 



In the strict adiabatic limit Plz = and we denote 
the respective matrix by U^ z ■ I n the first an d second 
crossings Plz might be different and accordingly we 
use the notations U^ z and U® z . The diagonal matrix 
ip — diag{<^+, if-} contains the dynamical phases that 
are accumulated in the upper and lower levels during 
the time between the two crossings, and we use the 
notation (p for the phase difference. The transposition 
operator T is required because in the second half of 
the cycle the roles of the lower and the upper states 
are interchanged. The expression for the total Landau- 
Zener transition probability is an interference of the 
two possible ways to get to the upper level, either by 
making the transition in the first crossing or in the 
second crossing: 



P 



(28) 



Physically p is the probability not to come back to the 
initial site. The strict adiabatic limit is p — 0. 

We turn now to calculate Q in leading order for 
Plz <C 1. The operator Q^z is obtained by integrating 
overX(t) = U{tflU{t) with U(t) = U° z (t). But in the 
second half of the cycle U(t) is given by Eq.(26), with 
U® z replaced by U^ z (t). Consequently we get 



l(t)dt 



Q2z 



[Te-^uSzf Qlz [Te- l *U? z 



The first term in this expression is Eq.(9) multiplied 
by Ao, with the p of Eq.(12), and with <f> = — 7r/2 cor- 



responding to the phase convention in Eq.(27). Then 
the second term in this expression becomes 



(1 — Sp) iSq 
—iSq +(1 — Sp) 



where 

Sp 
Sq 



+2P° Z + P° 




(29) 



(30) 
(31) 



It is easily verified that for A = 1 we have indeed 
an agreement with the restricted quantum-classical 
correspondence relation of Eq.(15) where p is given 
by Eq.(28). For A 7^ 1 the result for the expecta- 
tion value (Q) is in agreement with Eq.(25). For the 
variance we get 



Var(Q) 



'P L ° z + e^A 01 



(32) 



where Xq = A — 2A . This is a generalization of a re- 
sult that we had obtained using an adiabatic formalism 
in a previous publication [10]. The appearance of Aq 
instead of Ao reflects the definite site preparation at 
t = 0, and therefore it did not emerge in the adiabatic 
calculation where t = has no special meaning. Notice 
that the adiabatic approximation with Ao = Aq is for- 
mally obtained by replacing C/lz with U^%. 

One should realize that the interference expresses it- 
self differently in the expressions for p and for Var(Q). 
One may re-phrase this observation by saying that for 
a 3-site ring geometry, unlike the case of a 2-site geom- 
etry, there is no trivial relation between the counting 
statistics and the occupation statistics. 



7. Long time statistics of induced currents 

In this section we would like to consider the long time 
behavior of the counting statistics for either non-driven 
or periodically driven systems. In the latter context 
it is convenient to define U = U(T) as the one-period 
(Floque) evolution operator, and Q = Q(T) as the one- 
period counting operator, where T — 1 is the period. 
The interest is in the counting statistics after t periods. 
Accordingly Eq.(3) takes the form 
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Q(t)= f X{t')dt' = Q + U~ l QU + ■■■ +lf- t QU t 
Jo 

It should be clear that in this latter discrete version the 
operator Q (which is like flow per period), plays the 
same role as the operator X (flow per unit time). For 
this reason we are not going to duplicate the discussion 
below, and stick to continuous time notations. Having 
defined X as the current trough a specified bond (or 
more generally it would be replaced by the flow per- 
period), we can decompose it into a "DC" part and 
oscillatory part as follows: 

X(t) = ^\n)e i(En - Em)t I nm {m\ = X + SX(t) (33) 
where 



E n =E m 



n)I nm {m\ 



(34) 



In the absence of magnetic fields the stationary state 
of non-driven system is characterized by zero "DC" 
current and we get X = 0, unless the Hamiltonian has 
a degeneracy. For a periodically driven systems n are 
defined as the eigenstates of the Floque operator and 
in general we have 1^0. Accordingly 



Q(t)=tX+ J2 I") 



.1 - e 



i{E n -E m )t 



E n E m 



(m\ 



The non-zero elements of the oscillatory term are all 
off-diagonal, while those of X may be both diagonal and 
non-diagonal. However, without loss of generality we 
can choose the n basis such that X is diagonal. If the 
preparation is a superposition of Floque eigenstates, 
then both the average and the dispersion grow linearly. 
Accordingly, in general, the asymptotic value of the 
relative dispersion ^/Var(Q)/(Q) does not go to zero. 
In order to make it go to zero we have to especially 
prepare the system in an n eigenstate. For such special 
preparation {Q{t)} grows linearly while the dispersion 
is oscillating around a constant value: 



Var(Q) = J2 



2 \lnrn \ ' 
(E n — E m ) 



2 [1 - cos((E r -E m )t)] 



Var(Q) + oscillations 



(35) 



The time averaged value of the dispersion Var(Q) is 
given above by the first term (without the cosine). The 
time averaged value is not always of physical interest. 



If only two levels are involved, then Var(Q) drops to 
zero periodically. This applies also if several levels are 
involved, as long as their spacing differences are not 
resolved. The example in the following section demon- 
strates these considerations. 



8. Example for long time counting statistics 

The simplest way to illustrate the discussion of the 
previous section is to consider the counting statistics 
which is associated with the persistent current in a 
clean ring. For simplicity we consider N site ring. 
The position eigenstates are labeled as x = 0, 1, 2, 
while the momentum eigenstates are labeled as 
V = k n = (2n/N)n with n = 0, ±1, ±2, .... The Hamil- 
tonian is 



-cD - cD~ 



-2c cos (p) 



(36) 



where D — exp(— ip) is the one site displacement 
operator, and c is the hopping amplitude. The eigen- 
energies are E n = —2ccos(2-kii/N). Thanks to the 
n i — ► — n degeneracy of the spectrum the ring can sup- 
port a persistent current even in the absence of a mag- 
netic field. The velocity operator is a three diagonal 
matrix 



v = i[H, x] — icD — icD 1 = 2csin(p) 



(37) 



We measure the current through the 0^1 bond. Ac- 
cordingly 



i><oMo>(i|] 



X — —ic 
We realize that 

[X(t)] nm = ^cl[e lfc " -e-'*™] c 1 



(E n ~E m )t 



(38) 



(39) 



Upon time averaging only the n — m terms survive 
and they equal — i(l/N)2c sin(fc n ). Thus we have the 
identification 



(40) 



If the energies were equally spaced with some level spac- 
ing A = (2tt/N)vf, the motion of the particle would 
be strictly periodic. The period of the motion 2n/A 
would be the time to make one round along the ring. In 
such case one easily realizes that the variance Var(Q) 
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becomes zero at the end of each period. But if one 
takes the true dispersion into account, one realizes that 
this periodicity is not strict: after some time the quasi- 
periodic (rather than periodic) nature of the motion is 
exposed. The long time average of the fluctuating vari- 
ance can be calculated using Eq.(35), leading to: 



Var(Q)= Yl 



2\I„ 



^ „2 6 ( > 



Thus, even if a particle is prepared in a definite station- 
ary velocity eigenstate, still the counting at the end of a 
period does not yield a certain result. This uncertainty 
is related to the non-linearity of the spectrum. 



More generally we can define an "associated" phys- 
ical problem as follows. Assume that the particles are 
charged (for simplicity we set e = 1), and the current 
is driven by an electro-motive-force (EMF) which is in- 
duced by a vector potential 



A(x,t) = $(t)6{x-x ) 



(44) 



Then the rate of energy absorption is proportional to 
the current and to the EMF (-$) 



4x1 



(45) 



m m m 

dt dt <9$ 

Thus the transported charge Q implies energy absorp- 
tion 



Qabsorption = {H{t)-TL{0)) = -$ X Q 



(46) 



9. Counting, spreading, and dissipation 

The analysis of the long time behavior of the count- 
ing statistics along the lines of the previous section is 
not very illuminating once we turn to consider driven 
systems of greater complexity. Technically this is be- 
cause the quasi-energies, unlike E n become dense in 
the [0, 2n] interval. It is therefore more illuminating to 
observe that in the latter case the theory of counting 
is strongly related to the theory of spreading in real or 
in energy space. 

Let us assume that we have a ring, and that the 
current is measured through a section x. [In the tight 
binding model we can associate a location Xi with each 
site, and Eq.(38) is an example for an expression for 
the current through a section at xo < x < xi]. For a 
ring of length L it is convenient to re-define the current 
operator as: 



LJo 



dxX 



(42) 



which is essentially the velocity operator. Thus we get a 
relation between the counting operator and the spread- 
ing in real space 



Q(t) = !(*(«)- s(0)) 



(43) 



where x(t) is the non-periodic extension of the position 
operator. This procedure is nice but it should be clear 
that the re-definition of the current operator implies 
that possibly important information is lost. 



The energy absorption can be either positive or nega- 
tive, but if the system is chaotic it can be argued that 
its dispersion, and hence the average are growing as a 
function of time. The same applies to Q. Namely: 

Var(Q) =11 {I{t)l(t'))dtdt' (47) 
Jo Jo 

Thus both the counting statistics and the dissipation 
reflect the fluctuations of the current (Kubo). This 
point of view is quite powerful. Instead of thinking 
about "counting statistics" (which is quite abstract) 
we can think about "energy absorption" for which we 
have better intuition and better theoretical tools for 
analysis. 



10. Summary 

In closed geometries counting statistics does not 
obey the common stochastic point of view, but rather 
reflects the coherent nature of the quantum transport. 
For a single path crossing in a 2 site system the first 
and second moments coincide with the naive classi- 
cal expectation due to a restricted quantum-classical 
correspondence principle that can be established. 
But more generally, in multiple path geometries, the 
results are not as naively expected. 

In a double path geometry the particle has two paths 
of transport. The splitting ratio A would be in the range 
[0, 1] if the particle were classical. But in the case of co- 
herent transport A can be any number (either negative 
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or positive). We have explained what is the correct way 
to incorporate the splitting ratio A in the calculation 
of the counting statistics. 

During a double path crossing the probability ampli- 
tude of the particle is transported simultaneously via 
the two available paths. The coherent splitting is "ex- 
act" , and consequently Q has zero variance. If the par- 
ticle had finite probabilities to go either via the first 
path or via the second path, it would imply a non-zero 
variance (which is not the case). 

Possibly the most interesting result is the analysis 
of a full stirring cycle where the counting statistics 
becomes unrelated to the changes in the occupation 
statistics. In particular we showed that the interference 
of sequential Landau-Zener crossings is reflected differ- 
ently in the respective expressions for VarQ and for p. 

The RMS of the fluctuations of Q in a coherent stir- 
ring process grows like oc t and not like oc This 
linear increase can be avoided, or at least minimized, 
if there is control over the preparation of the system. 
Then we are left with a constant residual dispersion. 
For a driven system the counting statistics is related to 
the study of spreading and dissipation in energy space, 
and hence the growth of the variance constitutes a re- 
flection of linear response characteristics. 
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